Possible origin of the 0.5 plateau in the ballistic conductance of quantum point contacts 
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Abstract 

A non-equilibrium Green function formalism (NEGF) is used to study the conductance of a side-gated 
quantum point contact (QPC) in the presence of lateral spin-orbit coupling (LSOC). A small difference 
of bias voltage between the two side gates (SGs) leads to an inversion asymmetry in the LSOC between 
the opposite edges of the channel. In single electron modeling of transport, this triggers a spontaneous 
but insignificant spin polarization in the QPC. However, the spin polarization of the QPC is enhanced 
substantially when the effect of electron-electron interaction is included. The spin polarization is strong 
enough to result in the occurrence of a conductance plateau at 0.5 Go (Go= le'lh) in the absence of any 
external magnetic field. In our simulations of a model QPC device, the 0.5 plateau is found to be quite 
robust and survives up to a temperature of 40K. The spontaneous spin polarization and the resulting 
magnetization of the QPC can be reversed by flipping the polarity of the source to drain bias or the 
potential difference between the two SGs. These numerical simulations are in good agreement with 
recent experimental results for side-gated QPCs made from the low band gap semiconductor InAs. 



I. Introduction 



The signature of 1-D ballistic transport is the quantization of its conductance. The conductance 
exhibits plateaus at integral multiples of Go = 2e^/h [1, 2]. These plateaus are well understood in 
the framework of the one-elecfron theory [3]. In 1996, Thomas and co-workers observed an 
additional plateau at 0.7 Go in the AlGaAs/GaAs quantum point contact (QPC) in the absence of 
any magnetic field [4]. Since then, this anomalous plateau, usually referred to as the "0.7 structure", 
has been observed in both QPCs and relatively long ID quantum wires [5-9]. 

To understand its origin, the evolution of the 0.7 structure with temperature, applied magnetic field, and 
applied voltage has been studied experimentally. Typically, the 0.7 structure exhibits an anomalous 
temperature dependence showing increase of the conductance to Go as the temperature is lowered [5]. 
The zero-bias peak (ZBA) has been observed in the non-linear differential conductance and the 
appearance of the ZBA peak is related to the disappearance of the 0.7 structure [10]. In a parallel 
magnetic field, the 0.7 structure evolves smoothly into the Zeeman spin-split O.SGq plateau [4]. At the 
same time, the actual position of the 0.7 plateau may vary in the range [0.5 to 0.7] Gq depending on the 
electron density, the chaimel length and the lateral confining potential [11, 12]. 

The 0.7 structure has also generated intense theoretical efforts and several models have been proposed, 
but at present there seems to be no consensus on its origin. The Kondo model successfully predicts the 
temperature dependence and the ZBA for the 0.7 structure [13, 14], but it contradicts with the recent 
observation of the static spin polarization in a hole QPC [10]. Based on strictly one-dimensional 
Luttinger liquid state, a Wigner-cystal model has been used to explain the 0.7 structure and its 
temperature dependence [15]. However, a QPC is not necessarily a strictly ID system, so the Luttinger 
liquid theory may not be applicable here. Another model links the 0.7 structure to the spontaneous spin 
polarization in the QPC [16-19]. Since the static spin polarization has been experimentally found to be 
associated with the 0.7 structure [10], it deserves our attention. The local exchange energy can 
introduce a spin-spht energy gap. When Fermi energy lies in the spin-split gap, a static spin 
polarization gives a plateau in the range of [0.5 to 0.7] Gq. However, all these spin polarization models 
need to implant magnetic impurities or apply an initial small magnetic field [16-19] to trigger initial 
spin imbalance in the simulations. 

In a recent report [20], we presented experimental evidence of the 0.5 plateau in side-gated quantum 
QPCs, when the confining potential of the QPC is made sufficiently asymmetric. In these structures, 
the asymmetrical lateral spin orbit coupling (LSOC) is used to introduce the initial spin unbalance [20]. 
The QPCs were made from InAs (with a large infrinsic SOC, more than one order of magnitude larger 
than GaAs [21]) quantum-well (QW) structures with a two-dimensional electton gas (2DEG) in the 
well. Figure 1(a) is the SEM of a QPC device made on InAs QW structure. The QPC channel is created 
by negatively biasing the two side gates (SG) UG and LG. The confining potential of the QPC is 
confrolled by the bias voltages of the side gates and can be made asymmetric by applying unequal 
voltages to these gates. Figure 1(b) shows a representative plot of measured conductance at 4.2K of the 
InAs QPCs as a function of common-mode side gate voltage in the absence of any applied magnetic 
filed. An additional short plateau at conductance G ~ 0.5 Gq was observed when the confining potential 
was made highly aymmetic (AS). It was absent in the symmetric (S) case. Such a plateau is known to 



result due to Zeeman spin splitting of the lowest ID subband in an applied magnetic field. This happens 
when the Fermi level lies in the spin-split energy gap. The occurrence of a 0.5 plateau in the ballistic 
conductance is a signature of complete spin polarization. The observed 0.5 structure of Fig. 1(b) 
indicates spontaneous spin polarization of the conduction electrons. It was observed only when the 
transverse confining potential of the QPC was made asymmetric by appropriately adjusting the 
side-gate voltages. The 0.5 plateau was also observed when the confinement asymmetry was reversed 
by flipping the asymmetry of the gate voltages. 

In ref [20], a simple analytical model of the potential energy profile felt by the electrons in the central 
portion of the QPC was used in conjunction with a non-equilibrium Green function formalism (NEGF) 
to provide a theoretical explanation of the experimentally observed 0.5 plateau of Fig. 1(b). We showed 
that the 0.5 plateau is related to a spontaneous spin polarization induced by LSOC in a side-gated QPC. 
More precisely, the 0.5 Go conductance plateau appears in QPCs in the absence of any external 
magnetic field as a result of three ingredients: an asymmetric lateral confinement, a LSOC, and a 
strong electron-electron interaction. The asymmetry in the LSOC is required to trigger a small initial 
spin imbalance. This is true even in a single-electron description of carrier transport through the QPC 
[22]. However, the small spin imbalance is strongly enhanced when the effect of electron-electron 
interaction is taken into account to a point that near-perfect spin polarization in the QPC is reached 
resulting in the 0.5 plateau in the conductance. In this paper, we investigate in more detail the nature of 
the 0.5 plateau by considering more realistic potential energy profiles in side-gated QPCs using a 
self-consistent NEGF that includes the effect of space-charge in the structures. We study the effects of 
temperature and strength of the electron-electron interaction on the shape of the 0.5 conductance 
plateau. This paper is organized as follows. In the next section, we discuss in detail the origin of the 0.5 
plateau. In section III, we describe the NEGF approach used to calculate the charge densities of spin-up 
and spin-down electrons throughout the QPC. Section IV gives the results of our numerical simulations. 
Finally, Section V contains our conclusions. 



II. The Origin of the 0.5 Plateau 

We model the QPC using the configuration shown in Fig. 2 where the white region represents the mesa 
etched quantum wire with openings on either side. The gray areas represent the etched isolation 
trenches that define the dimensions of the QPC. There are four contacts connected to the QPC device, 
source, drain and two SGs. Symmetric and asymmetric SG voltages can be applied. Since the QPC of 
Fig. 1(a) is made from a nominally symmetric InAs QW, spatial inversion asymmetry was assumed to 
be negligible along the growth axis (z axis) of the QW and the corresponding Rashba spin-orbit 
interaction was neglected. The Dresselhaus spin-orbit interaction due to the bulk inversion asymmetry 
in the direction of current flow is also neglected. The only spin-orbit interaction considered is the 
LSOC due to the lateral conflnement of the QPC channel, provided by the isolation trenches and the 
bias voltages of the side gates. To understand the origin of the 0.5 plateau, we begin with the 
single-particle Hamiltonian, which is given by 



where i/^ = — -ip^ + p^j+U{x,y) , P is the intrinsic SOC parameter, a is the vector of Pauli spin 

matrices, and fi^^ is the effective magnetic field, which is induced by the LSOC. The 2DEG is 

assumed to be located in the (x, y) plane, x being the direction of current flow from source to drain and 
y the direction of confinement of the channel. t/(x, y) is the confinement potential, which includes the 
potential inttoduced by contact gates and conduction band discontinuity b£^{y) at the InAs/air 

interface. 

Figure 3(a) gives a schematic representation of the confining potential along y direction when a 
symmetric side-gate voltage is apphed. The effective magnetic field has exactly the same 

magnitude but opposite directions at the opposite transverse edges of the QPC. Moving electrons with 
opposite spins experience opposite SOC forces that leads to an accumulation of opposite spins at the 
opposite transverse edges. The spin-up is the majority spin species on edge I of the QPC and the 
minority species on edge II. The difference of spin density is anti-symmefric about j = wi/2 giving zero 
net polarization. 

When an asymmetric SG voltage is applied on the QPC, the potential profile changes from the 
symmetric dashed line to the asymmetric full hne as shown in Fig. 3(b). The spin-up population on the 
left edge I exceeds the spin-down one on the right edge II. This results in a net spin-up polarization, 
which is the initial imbalance between spin-up and spin-down electrons induced by the asymmetric 
LSOC. It can be shown that when the asymmefry between Vsgi and Vsg2 is reversed, the direction of 
spin polarization is reversed. The strong repulsive Coulomb e-e interaction enhances this imbalance. 
As a result, the spontaneous spin polarization can reach nearly 100% in the regime of single-mode 
ttansport and the 0.5 conductance plateau can result. 

III. The NEGF Model for Side-Gated QPC 

The conductance through the QPC was calculated using a NEGF method under the assumption of 
balhstic transport [23]. Within the NEGF, the Green function associated to the QPC is then a 

( 2N^N^ X 2N^Ny ) matrix (where A^x and A'y are the number of grid points along and perpendicular to 

the direction of current flow and the factor two is needed to distinguish the degree of spin polarization) 
and is given by 

G(i?) = (i?/-//-S,-S„-S^.)->, (2) 

where 2^ and are the self-energy terms representing the coupling of the source and drain 
contacts [23] and the electton-electron interaction self-energy. We used a Hartree-Fock 



approximation following Lassl et al. [19] to include the effects of Coulomb e-e interaction in the QPC. 
In this approach, the interaction between two electrons located at (x,>')and (x',y') is modeled using 

the contact potential 



(3) 



where y is the interaction strength equal to a few times of I Im . 



The interaction self-energy with spin o is then given by 



Ktix,y) = 7n_^ix,y), 



(4) 



where n_^{x,y) is the density of electrons with spin -a. 

The interaction self energy 'Z^^(x,y) is different for the two spin populations injected from the 

contacts. A spin-up electron encounters a potential which is proportional to the density of spin-down 
electrons, and vice versa. This leads to a repulsive interaction between electrons with opposite spin 
directions. Any external source leading to an imbalance between the density of spin-up and spin-down 

electrons is increased by the addition of the self-energy termS^,(x,>') . In our case, it is the asymmetric 
LSOC, which leads to the initial imbalance. 

Once H, Xs> Zd and Zint are known, the Green function (G) can be calculated from Eq. (2) and all the 
other quantities of interest can be found out from the following set of equations. 

1) The broadening mattices for the source and drain contacts 



r, (E) = io:, -lis'), (E) = iiXn - Sz, 



(5) 



2) The spectral functions of the source and drain contacts 



A,iE) = Gr,{E)G\ A^(E) = Gr^{E)G\ 



(6) 



3) The density matrix 



00 



-dE, 



(7) 



-co 



where^ and fo are Fermi-Dirac distributions in the source and drain contacts, respectively. 



5) The transmission coefficient 



T{E) = Trace[TfiTj,G% 



(8) 



6) The current through the QPC under the assumption of ballistic transport 



I = UnE)\f,{E)-f^{E)\dE. (9) 

— QO 

Based on Eq.(5), Eq.(6) and Eq.(7), we can get the density matrix p from the Green ftinction mattix G. 
The diagonal elements of the density matrix determine the spin-up density and spin-down density «; 
and hence the interaction self-energy Xint from Eq.(4). Since the interaction self-energy depends on the 
Green function matrix and vice-versa, an iterative procedure is required to obtain the final results. 



IV. Numerical Results 



We consider a side-gated QPC made from InAs QW structure with a 2DEG in the well. The low band 
gap semiconductor InAs has a large intrinsic SOC. The effective mass in the InAs channel was set 
equal to m*= 0.023/«o, where wq is the free electron mass. Unless otherwise stated, all calculations were 
performed at a temperature T= 4.2 K. Following Lassl et al. [19], the strength of the parameter y in the 

interaction self energy was set equal to3. 7 h^/lm* . The strength of the parameter in the LSOC was 

set equal 200 [20]. Unless otherwise stated, the geometrical parameters /i, k, wi and W2 were 
selected to be 68, 36, 48 and 16 nm, respectively. These parameters are smaller than the experimental 
values of the QPC shown in Fig.l (a) and were chosen to reduce computational time. The potential at 
the source was set equal to OV and the one at the drain to 0.3 mV in all simulations. An asjmimetry 
in the potential of the SGs was introduced by taking Fsgi = 0.2 V + Fsweep and = -0.2 V + Fsweep^nd 
the conductance of the constriction was studied as a function of the sweeping (or common mode) 
potential, Fsweep- The Fermi energy was equal to 106.3 meV in the source contact and 106meV in the 
drain contact, ensuring single-mode transport through the QPC. 



At the interface between the rectangular region of size W2X/2 and vacuum, the conduction band 
discontinuity was modeled as 



at the bottom interface, and 



AE. 



1 + cos— y 



l-HCOS- 



(10) 



(11) 



at the top interface to achieve a smooth conductance band change, where d was selected to be in the nm 
range to represent a gradual variation of the conduction band profile from the inside of the quantum 
wire to the vacuum region. A similar grading was also used along the walls going from the wider 
portion of the channel to the central portion of the QPC, as shown in Fig. 2. The gradual change in 
AE Xy) is responsible for the LSOC triggering the spin polarization of the QPC in the presence of an 

asymmetry between V^gi and Vsg2, as discussed below. The parameter d appearing in Eqns. (10) and (11) 
was set equal to 1.6 nm. Similar results were obtained when AE^(y) was linearly changed at the 

interface and d was set equal to 0.8 nm and 1 .2 nm. The conductance of the QPC was then calculated 



using the NEGF outlined in section II using a non-uniform grid configuration containing more grid 
points at the interface of the QPC with vacuum [24]. 

Figure 4 is our main result which shows plots of the conductance of the QPC as a function of Fs^eep for 
symmetric (AFsg= Vsgi-Kgi^ 0) and asymmetric (AFsq = Fsgi-Fsg2= 0.4 V) confinements. The 
dashed curve is the conductance calculated with the symmetric confinement and only one plateau at 
2e^/h is observed. The oscillation in the conductance for Fsweep > is a result of multiple reflections 
between the ends of the central rectangular portion of the QPC. The fiill curve labeled "G-^ + G;" is the 
conductance with asymmetric confinement and clearly indicates the presence of a plateau in 

conductance around e^/h besides the normal 2e^/h plateau. Not shown here, the contributions of 

and G| to the conductance was found to be switched when the polarity of A Fsg was flipped. Figure 

4 also shows the contribution of the majority and minority spin bands as a function of Fjweep- G^ and G| 
are found to be decreasing and increasing, respectively, near Fsweep = V leading to a negative 
differential region in G| + G|. This feature is quite common in the numerical simulations of QPCs 
[16,17,19,25]. It is partly due to the effect of multiple reflections at the edges of the narrow portion of 
the QPC but also depends very strongly on the exact shape of the potential energy landscape in the 
QPC [16]. As shown in Fig. 5, the size of the conductance modulation A on the plateau gets smaller 
when the aspect ratio W2//2 is closer to unity, in agreement with the experimental results shown in Fig. 1 
for which W2//2 ~ 1 • 

Figure 6 shows that the 0.5 plateau is rather sensitive to the choice of the parameter y, but is otherwise 
robust. If the constriction of the QPC is narrow enough for transport to be single-mode, the 0.5 plateau 
is fairly robust as a function of temperature, as shown in Fig.7. Even at a temperature of 40 K, a 
remnant of the 0.5 plateau can still be seen. The difference in conductance G^and G| in the transition 
region from the cut off region and the first quantized plateau 2e^/h can be understood by plotting the 
potential energy profile for the majority and minority spin bands. Figure 8 shows a plot of potential 
energy profile and in the central portion of the QPC for Fsweep = V, where U^{x, y) = U{x, y) + 
Y,mt and Ui(x, y) = U(x, y) + Y.mt- The potential energy profile {/; has a camelback shape which 
prevents the flow of minority spins through the QPC, whereas U-^ has a saddle shape. The Fermi level is 
then in between the local maximum (point A in Fig. 8(a)) and local minimum (point B in Fig. 8(b)) felt 
by the minority and majority band, respectively. The majority spin band is propagating through the 
channel but the minority spin band is evanescent. The difference in potential energy between points A 
and B located at the center of the rectangular portion of the QPC is equal to 0.033eV which is about 
100 times = 0.36 meV for T = 4.2 K. This explains why the 0.5 plateau can still be observed at 
temperature as high as 40K in Fig.7. With further increase in Fsweep, the Fermi level is now above both 
maxima of the potential energy profile f/| and {7| in the central portion of the QPC, both channels are 
propagating and the difference between f/| and f/| disappears gradually. The conductance then reaches 
a maximum of 2e^/h. 

A plot of «| (x,y) - «| (x,y) for Vswccp = is shown in Fig. 9 to illustrate that the spin imbalance is the 
strongest in the central (rectangular) portion of the QPC. To better illustrate the efficiency of spin 
injection through the QPC, we have plotted in Fig. 10 the spin polarization 



a(x,y)={n^(x,y)-n^(x,y))/{n^(x,y) + n^(x,y))- Its maximum value was found to be 81% but it can 
be tuned closer to 100 % by changing the QPC structural parameters. For comparison, we plot in 
Figures 1 1(a) and 1 1(b) the quantities (x,y) - (x,y) and a(x,y) when the parameter y is set equal to 
zero, i.e., the effects of electron-electron interaction are neglected. In this single electron picture of 
carrier transport, the imbalance induced by the asymmetry in LSOC triggers an accumulation of 
opposite species on the opposite edges of the central portion of the QPC [22]. This spin imbalance 
leads to an asymmetry in a(x,y) along the width of the chaimel, as shown in Fig.l 1(b). In this figure, 
the maxknum of a(x,y) is about a factor 400 below the maximum observed in Fig. 10. This stresses the 
importance of electron-electron interaction in the chaimel due to the substantial difference between the 
potential energy profiles f/| and {7| acting on the two different spin species in the central portion of the 
QPC. 

In summary, NEGF simulations of spontaneous spin polarization in side-gated QPCs present the 
following picture. Below threshold, the Fermi level in the source is below the minimum of the potential 
energy profile for either spin band and the conductance is close to zero. When the confining potential 
of the QPC is symmetric, opposite spin accumulations induced by LSOC at the transverse edges cancel 
each other and there is no net spin polarization. When the confinement potential is made asymmetric by 
applying a small potential difference between the two side gates, the asymmetric LSOC triggers a small 
imbalance between the majority and minority spin bands. When this imbalance is fed back into the 
local e-e interaction self-energy term, the difference in the potential energy profiles felt by the two 
different spin chaimels becomes quite drastic [26]. In addition, we have shown that the shape of the 
anomalous plateau is quite robust but otherwise sensitive to the value of the parameter y in Eq.(4), i.e., 
the exact form used to model the effect of exchange-correction potential energy. Extensive conductance 
measurements aroimd the 0.5 plateau of QPC could therefore be used to refine the details of the 
theoretical treatment and understand better the importance of e-e interaction. 

V. Conclusion 

The use of asymmetrically biased lateral QPCs offers an all-electrical way to generate highly 
spin-polarized current avoiding the need for ferromagnetic contacts or external magnetic field. The use 
of two QPCs in series with a submicron long channel in between whose width could be tuned by two 
additional side gates could therefore pave the way for the first demonstration of an all-electrical 
realization of the Datta-Das SpinFET. The same approach could provide an all electrical means to 
realize semiconducting quantum computing gates. 
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Figure Captions 



Figure 1(a) Scanning electron micrograph (SEM) of a side-gated InAs QPC. SGI and SG2 are side 
gates. Dark areas are trenches cut by wet etching. Light areas are the wafer surface with the 2DEG 
underneath, (b) Conductance of InAs QPC at 4.2K in the absence of applied magnetic filed. The 0.5 
conductance plateau is observed when the confining potential is asymmetric (AS). It is absent when the 
confining potential is symmetric (S) [20]. 

Figure 2: Schematic of the QPC configuration used in the numerical simulations. The geometrical 
dimensions h, h, wi and W2 are selected to be 68, 36, 48 and 16nm, respectively. In all simulations, Vs = 
OV, Fd = 0.3mV. 

Figure 3: (a) Schematic representation of the confining potential energy along y direction at x = /i/2 
when symmetric SG voltages are applied; (b) Same for asymmetric SG voltages. 

Figure 4: Conductance of a QPC as a function of Fsweep with the geometrical dimensions shown in Fig.2. 
The SG biasing parameters for the symmetric confinement case are Fsgi = Vsgi = Vsweep', for the 
asymmetric case, Fsgi = 0.2V + Fsweep and Vg^ = -0.2V + Fsweep- The temperature is set equal to 4.2K. 
The top dashed curve corresponds to the symmetric confinement. The individual contributions of the 
two spin bands to the full curve "G^ + G" are is shown when asymmetric confinement is applied in the 
simulations. The 0.5 plateau is clearly visible then. 

Figure 5: Conductance of a QPC as a fiinction of Fsweep close to the 0.5 plateau for different choices of 
the aspect ratio W2//2. The parameter y is set equal to 3.7 in units of I Im . = OV, = 0.3mV, Vsgi 
= 0.2V + Fsweep, Vsgi = -0.2V + Fsweep- The temperature is 4.2K. The geometrical dimensions h, w\, and 
W2 are selected to be 68, 48 and 16nm, respectively. From left to right, the curves correspond to h equal 
to 36, 32, 26, and 24 nm, respectively. The corresponding value of the conductance modulation A on 
the plateau is equal to 0.29, 0.24, 0.17, and 0.12 (lelK), respectively. 

Figure 6: Dependence of the 0.5 conductance plateau on parameter y (in units of ti^llrn). The 
biasing parameters of the QPC are = OV, Vi = 0.3mV, Fsgi = 0.2V + Fsweep, V^^ = -0.2V + Fsweep- The 
temperature is 4.2K. The geometrical dimensions h, h, wi and W2 are selected to be 68, 36, 48 and 
16nm, respectively. 

Figure 7: Temperature variation of the 0.5 conductance plateau for a QPC with the geometrical 
configuration shown in Fig.2. The biasing parameter of the QPC are Vs = OV, = 0.3mV, V^gi = 0.2 V + 
^sweep; ^sg2 = -0.2V + Vsv/eep- The parameter y is set equal to 3.7 in units of 1 2m . The geometrical 
dimensions 1\, I2, W\ and W2 are selected to be 68, 36, 48 and 16nm, respectively. 



Figure 8: (a) Potential energy profile for minority spin band in the central portion of QPC. (b) Same for 
majority spin band. The parameter y is set equal to 3.7 in units of ti^ I2rn , Fswccp ^ OV, and the 
temperature T = 4.2K. The geometrical dimensions l\, h, w\ and W2 are selected to be 68, 36, 48 and 
16nm, respectively. There is a difference of 33.6meV between the points labeled A and B. 

Figure 9: Spatial variation of the difference of the density for majority and minority spins across the 
QPC for V, = OV, Vi = 0.3 mV, Fsgi = 0.2V, V,^ = -0.2V (i.e., F^weep = OV). The temperature 7= 4.2K 
and y = 3.7 in units of llm . The geometrical dimensions h, h, wi and W2 are selected to be 68, 36, 
48 and 16nm, respectively. 

Figure 10: Spatial dependence of the spin polarization a(x,y) = {n^(x, y) -n^(x, y))/{n^(x,y) + n^(x, y)) 
in the central portion of the QPC for V, = OV, = 0.3mV, Fsgi = 0.2V, V^gi = -0.2V (i.e., Fsweep = OV). 
The temperature T = 4.2K and y = 3.7 in units of I Im . The geometrical dimensions h, h, wi and 
W2 are selected to be 68, 36, 48 and 16nm, respectively. 

Figure 11: Plot of (a) the spatial variation of the difference of the density for majority and minority 
spins across the QPC and (b) the spatial dependence of the spin polarization 
a(x, y) = [n^ (x, y) - (x, y))l{n^ (x, y) + (x, y)) in the central portion of the QPC when the effects of 

the electron-electron interaction in the QPC are neglected (y = 0). The biasing conditions, physical and 
device parameters are the same as in Figures 9 and 10. 
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